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We analyze thermonuclear and pycnonuclear fusion reactions in dense matter containing atomic 
nuclei of different types. We extend a phenomenological expression for the reaction rate, proposed 
recently by Gasques et al. [rO ] for the one-component plasma of nuclei, to the multi-component 
plasma. The expression contains several fit parameters which we adjust to reproduce the best 
microscopic calculations available in the literature. Furthermore, we show that pycnonuclear burning 
is drastically affected by an (unknown) structure of the multi-component matter (a regular lattice, 
a uniform mix, etc.). We apply the results to study nuclear burning in a 12 C 16 mixture. In this 
context we present new calculations of the astrophysical ^-factors for carbon-oxygen and oxygen- 
oxygen fusion reactions. We show that the presence of a CO lattice can strongly suppress carbon 
ignition in white dwarf cores and neutron star crusts at densities p > 3x I0 9 g cm -3 and temperatures 
T < 10 s K. 



I. INTRODUCTION 

Nuclear reactions are most important for the physics 
of stars. They determine hydrogen burning in main- 
sequence stars, helium burning in red giants, and car- 
bon, neon, and oxygen burning at later stages. They 
determine also nucleosynthesis in shock driven stellar ex- 
plosions, such as type II supernovae, as well as ignition 
and burning in accreting stars. 

Steady-state and explosive thermonuclear carbon 
burning during late stellar evolution pj and in shock 
fronts of type II supernovae Q is governed by the 
12 C+ 12 C and possibly by the 12 C+ 16 fusion processes. 
Similarly, thermonuclear oxygen burning is mainly de- 
termined by the 16 0+ 16 and possibly by the 16 O+ 20 Ne 
reactions The ignition and nucleosynthesis during 
these burning phases critically depend on the initial fuel 
abundance and on the thermonuclear reaction rates. 

In high-density cores of white dwarfs and crusts of neu- 
tron stars, the thermonuclear reactions are strongly af- 
fected by plasma effects (especially important for carbon 
ignition in cores of accreting massive white dwarfs for 
triggering type la supernova explosions). The ignition 
conditions are sensitive to the 12 C and 16 O abundance 
and to the fusion reaction rates Q|. 

Carbon ignition has also been suggested as a trigger of 
superbursts in surface layers of accreting neutron stars 
However, the required ignition conditions seem to 
disagree with the observed superburst lightcurves @ . Al- 
ternative explanations are presently being discussed, such 
as carbon ignition in the crust of an accreting strange 
(quark) star, to accommodate the observed light curve 
characteristics Q. While in most of the scenarios pure 



carbon burning dominates the energy production, igni- 
tion conditions and associated nucleosynthesis are af- 
fected by the presence of other elements. 

Pycnonuclear burning occurs in dense and cold cores of 
white dwarfs and in crusts of accreting neutron stars 
0, ^3 • Theoretical formalism has been mostly developed 
for pycnonuclear reactions between equal nuclei, but one 
often needs to consider a multi-component matter, for 
instance, carbon-oxygen cores of white dwarfs. 

In a previous publication |l l| we have focused on fu- 
sion reactions between equal nuclei in a one-component 
plasma (OCP) of atomic nuclei (ions). In the present 
work we expand the study towards a multi-component 
plasma (MCP). The problem has two aspects; the first 
one is associated with the underlying nuclear physics 
while the second one is concerned with the plasma 
physics. The nuclear part deals with the reliable de- 
termination of astrophysical S-factors at stellar ener- 
gies. These energies are low (typically, lower than a 
few MeV), in particular if compared to the presently ac- 
cessible range of low-energy fusion experiments. This 
prevents direct measurements of S'-factors at laboratory 
conditions. Thus, one needs to calculate the S'-factors 
theoretically and use these results to extrapolate mea- 
sured S-factors towards lower stellar energy range. In 
Section|n]we present calculations of the S-factors for the 
two reactions of astrophysical importance, 12 C+ 16 and 

16Q + 16Q 

The plasma physics problem consists in calculating the 
Coulomb barrier penetration in nuclear reactions taken 
into account Coulomb fields of surrounding plasma par- 
ticles. These fields modify the reaction rates and lead 
to five nuclear burning regimes [jj (two thermonuclear 
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regimes, with weak and strong plasma screening; two py- 
cnonuclear regimes, for zero-temperature and thermally 
enhanced burning; and the intermediate regime). These 
regimes are described in Section IlIII their validity condi- 
tions are specified in Section flll Al In ordinary stars, nu- 
clear burning proceeds in the weak screening thermonu- 
clear regime [3 113 • The foundation of the theory of 
thermonuclear burning with strong plasma screening was 
laid by Salpeter The strict theory of pycnonuclear 
burning was developed by Salpeter and Van Horn @- 
References to other works can be found in Section IIIII 
and in [Tl| . 

In Section ITTT1 we analyze calculations of Coulomb bar- 
rier penetration in MCP and propose a phenomenologi- 
cal expression for a reaction rate valid in all five regimes 
for any non-resonant fusion reaction. This expression 
accurately reproduces well known results in thermonu- 
clear regimes and gives a reasonable description of the 
Coulomb tunneling problem in other regimes. It is im- 
portant for incorporating the plasma physics effects into 
computer codes which simulate nucleosynthesis, espe- 
cially at high densities in compact stars, such as white 
dwarfs and neutron stars (see above). In Section llVI we 
illustrate the results of Sections HTl and ITTT1 bv analyzing 
nuclear burning in 12 C 16 mixtures, with the emphasis 
on the carbon ignition curve. 



II. ASTROPHYSICAL S-FACTORS FOR 
CARBON-OXYGEN MIXTURES 

In order to study nuclear burning in dense stellar 
carbon-oxygen matter fSection IIVI) we need fusion cross 
sections (or associated astrophysical factors) for three re- 
actions, 12 C+ 12 C, 12 C+ 16 0, and 16 0+ 16 0. For calcu- 
lating the cross sections we employ the one-dimensional 
barrier penetration (BP ) formalism [l5l ] and adopt the 
Sao Paulo potential [IE 111 13 13 

to describe the real 
part of the nuclear interaction Vn {r, E) : 

V N (r, E) = V SP (r, E) = V F (r) e" 4 ^ . (1) 

Here, Vf(t) is the density-dependent double- folding po- 
tential, c is the speed of light, E is the particle collision 
energy (in the center-of-mass reference frame), v is the 
local relative velocity of two nuclei 1 and 2, 

v 2 (r, E) = {2/n)[E- V c (r) - V N (r, E)] , (2) 

Vc(r) is the Coulomb potential, and /j, is the reduced 
mass. 

In this paper, we adopt the two-parameter Fermi (2pF) 
distribution to describe the nuclear densities. The radii of 
these distributions are well approximated by the formula 
Rq = 1.31A 1 / 3 - 0.84 fm The 12 C and 16 diffuse- 
ness was taken to be 0.56 fm and 0.58 fm, respectively. 
These values were extracted from heavy-ion elastic scat- 
tering data at sub-barrier and intermediate energies, by 



applying an unfolding method involving the Sao Paulo 
potential (see Refs. pol I21T |22| for details). 

Usually, fusion cross sections a(E) at low energies, typ- 
ical for astrophysical conditions, are expressed in terms 
of the astrophysical S-factor 

S(E) = a{E) E c 27r " , (3) 

where r\ = [ZxZ^e 1 /K) yjij (2E) is the familiar Gamow 
parameter; Z\ and Zi are charge numbers of the nuclei. 
This parameterization removes from the fusion cross sec- 
tion the strong non-nuclear energy dependence 0, |2^| 
associated with Coulomb barrier penetration. If S(E) is 
a slowly varying function of E, it can be extrapolated to 
lower energies relevant to stellar burning. 

The S-factors for all three reactions versus E are shown 
in upper panels of Fig. ^ Solid lines are theoretical cal- 
culations using the BP model, while symbols are experi- 
mental data. The results for the C+C reaction have al- 
ready been discussed in Ref. 11] and are presented here 
for completeness of consideration. The data for this re- 
action are taken from Patterson et al. [24j|. Mazarakis 
and Stephens (2{|, Hig h and Cujec [2|J, Rosales et al. 
H3, Kettner et al. Hf, and Becker et al. The data 
for the C+O reaction are taken from Christenscn et al. 
[3fj| and Cujec and Barnes [3lJ . Finally, the data for the 
O+O reaction are from Hulke et al. l33|. Kuronen et al. 
[33j], Wu and Barnes [34|], Kovar et al. |35|. and Thomas 
et al. [3(|. Lower panels in Fig.^show the Gamow-peak 
energy ranges versus the temperature of stellar matter in 
the thermonuclear burning regime; they will be discussed 
in Section IhTaI 

The sub-Coulomb-barrier resonances exhibited in the 
C+C data at E < 6 MeV and in the C+O data at E < 7.7 
MeV cannot be reproduced in the framework of the BP 
model. However, the model provides an average descrip- 
tion of the light- and heavy-ion fusion at energies below 
and above the barrier. Depending on the nuclear poten- 
tial it also gives a satisfactory parameter-free description 
of the energy dependence of the S'-factor, which seems 
reasonably accurate for extrapolating the experimental 
data into the stellar energy range. 

The data sets presented by the different groups are not 
in perfect agreement. For instance, the two C+O data 
sets agree in average, but disagree within a factor of ~ 2 
for lowest E. The overall agreement between the theory 
and the data is ~50%. However, in the low-energy re- 
gion, the slope of the measured cross section reported in 
Ref. [3(j does not follow the calculated S'-factors. It is 
difficult to predict where the data will lie in the energy 
range E < 4 MeV. Regarding the O+O reaction, the dis- 
crepancies between the different experimental results at 
sub-barrier energies are around a factor of 3. The overall 
agreement between the data and the theory is ~30%. At 
the lowest measured energy the difference between the 
data and the theory are at most a factor of 3. Further 
experiments at lower energies would help in verifying the 
validity of the predicted fusion cross sections. Never- 
theless, it is important to highlight that the BP model 



3 




FIG. 1: (color online) Top: Astrophysical S-factors as a function of the center-of-mass energy E for the 12 C+ 12 C, 12 C+ le O, 
and le O+ 16 reactions. The solid lines correspond to the BP model calculations while the various symbols are experimental 
data. Bottom: Gamow-peak energy ranges for these reactions in the thermonuclear regime versus the temperature T of stellar 
matter. See text for details. 



does not contain any free parameter. In this sense, the 
5-factor calculations do not represent a fit to the experi- 
mental data and can be considered as a useful tool to pre- 
dict average non-resonant low-energy cross sections for a 
wide range of fusion reactions. For many astrophysical 
reactions, such a description gives a reasonable estimate 
because the formalism of stellar reaction rates relies on 
the S'-factor averaged over an entire Gamow-peak range. 

The values of S(E) calculated up to E < 20 MeV can 
be fitted by an analytic expression 



TABLE I: The coefficients A\,. . . ,As and Eq in the fits ex- 
pression (3) f° r the S-factors of the C+C, C+O and O+O 
reactions. 



Reaction 


Eo 


Ai 


A 2 


A 3 


A 4 


A 5 


12 C+ 12 C 

12 C+ 16 
16q + 16 


6.946 
8.290 
10.52 


34.75 
44.32 
56.16 


-0.552 
-0.561 
-0.571 


-2.131 
-1.480 
-1.160 


-0.625 
-0.910 
-1.044 


0.0315 
0.0387 
0.0366 



S(E) = exp(A 1 +A 2 AE 

A 3 + A A AE 

l + e- AE 



+ 



A A_E 2 \ 

^ J MeV barn, 



(4) 



III. NUCLEAR FUSION RATE 



A. Physical conditions and reaction regimes 



where AE = E — Eo; the center-of-mass energy E and 
the fit parameter E are expressed in megaelectron- volts. 
Table|Ugives the fit parameters A\, . . . , A5 and Eq for the 
C+C, C+O and O+O reactions. The maximum formal 
fit errors are 7.2% at E = 19.8 Mev for C+C; 6.3% at 
E = 7.5 MeV for C+O; and 3.9% at E = 8.2 MeV for 
the O+O reaction. The 5-factor for the C+C reaction 
was fitted in a previous paper by a slightly different 
expression with approximately the same accuracy. We 
have fitted the same data by the new expression Q) for 
completeness of consideration. The two fits are nearly 
equivalent. 



Let us consider a stellar matter which consists of ions 
and electrons. We assume that the ions are fully ion- 
ized and the electrons form a uniform electron back- 
ground. We study a multi-component mixture of ion 
species j = 1, 2, . . ., with atomic numbers Aj and charge 
numbers Zj. Let rij be the number density of ions j. The 
total number density of ions is n = JV rij ; the electron 
number density is n e = JV Zjrij. For an OCP of ions 
the index j will be omitted. The number density rij can 
be expressed through the mass density p of the matter as 
rij = Xjp/(Ajm u ), where Xj is the mass fraction of ions 
j, and m u = 1.66 5 5 x 10 - 24 g is the atomic mass unit. In 
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FIG. 2: (color online) Temperature-density diagram for a 
12q16q ma tt er . Straight lines show the temperature Ti of 
the appearance of ion liquid, the melting temperature T m of 
ion crystal, and the ion plasma temperature T v . Solid lines 
refer to pure carbon matter, long dashes to CO matter with 
equal particle fractions of C and O, and short dashes to pure 
oxygen matter. Three shaded strips show the regions impor- 
tant for C+C burning (in pure carbon matter), C+O burning 
(in CO mixture) and O+O burning (in pure oxygen matter). 
Each strip is restricted by upper and lower lines along which 
the burning time equals 1 year and 10 6 years, respectively 
(see the text for details). 



a not too dense matter the total mass fraction contained 
in the nuclei is Xm = J2j Xj ~ 1. At densities higher 
than the neutron drip density (~ 4 x 10 11 g cm~ 3 ), the 
matter contains also free neutrons; the total mass frac- 
tion contained in the nuclei is then Xn < 1. It is also 



useful to introduce the fractional number Xj 



nuclei j among other nuclei, with ^ 



rij/n of 



i(Z), 



P = 



i u n(A) 



1. Generally, 

X j/ A 3 



Xn ^ZiX-ijAi 

(Z) = £>i% (A) = J2x j A j , (5) 
i 3 

where (Z) and {A) are the mean charge and mass number 
of ions, respectively. 

Let us also introduce the Coulomb coupling parameter 
for ions j, 



Zje 2 _ z) 1 3 e 2 



(6) 



1/3 



, aj — Aj a e , 



where T is the temperature, fee is the Boltzmann con- 
stant, a e is the electron-sphere radius, and a.j is the ion- 
sphere radius (a radius of a sphere around a given ion, 
where the electron charge compensates the ion charge). 
Therefore, Tj is the ratio of a typical electrostatic en- 
ergy of the ion to the thermal energy. If Tj <C 1 then 
the ions constitute an almost ideal Boltzmann gas, while 
for Tj > 1 they are strongly coupled by Coulomb forces 
(constitute either Coulomb liquid or solid). The trans- 
formation from the gas to the liquid at Tj ~ 1 is smooth, 
without any phase transition. The solidification is real- 
ized as a weak second-order phase transition. According 
to highly accurate Monte Carlo calculations, a classical 
OCP of ions solidifies at T » 175 (see, e.g., Ref. [3fl|V 

It is useful to introduce the mean ion coupling pa- 
rameter (r) = ^2jXjTj. The plasma can be treated as 
strongly coupled if (T) > 1. This happens at T < Tj, 
where 



k B Ti = J2(Z]e 2 / aj ) Xl = k B T{T). 



(7) 



For low temperatures T <C T p <C Ti ion motion can- 
not be any longer considered as classical but should be 
quantized. Here, T p is the Debye (plasma) temperature 
associated with the ion plasma frequency oj p (a typical 
frequency of ion vibrations in Coulomb crystals - see, 
e.g., Ref. i), 



hhJp 

fc R 



E 

3 



AjWL u 



(8) 



The most difficult problem of a strongly coupled MCP 
at low temperatures consists in understanding its actual 
state. Extensive Monte Carlo simulations [37j of the 
freezing of a classical OCP indicate that it can freeze into 
imperfect body-centered cubic (bcc) or faced-centered cu- 
bic (fee) microcrystal (or microcrystals). Unfortunately, 
publications on reliable simulations of freezing of an MCP 
are almost absent. Evidently, the cold MCP is much more 
rich in physics than the OCP. It can be an MCP regular 
lattice or microcrystals (with defects); or an amorphous, 
uniformly mixed structure; or a lattice of one phase with 
random admixture of other ions. One cannot exclude an 
ensemble of phase separated domains. For the sake of 
completeness, we will consider different possibilities. An 
MCP obeys the linear mixing rule with the high accu- 
racy. Accordingly, the difference in energies of the in- 
dicated states is very small and is a subject of vigorous 
debates (e.g., [3^| and references therein). It is possible 
that a low-temperature MCP can be in different states 
depending on the history of its formation in a star with 
decreasing temperature. 

To make our consideration less abstract, we will apply 
it to a 12 C 16 mixture. The appropriate temperature- 
density diagram is shown in Fig. We present the 
temperatures 7], T m and T p for a pure carbon matter 
(xq = 1 , solid lines) , for a mixture of equal amounts of C 
and O nuclei (xc = 4, long dashes) and for a pure oxygen 
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matter (xc = 0, short dashes). The melting temperature 
of the CO mixture is taken to be T m = TJ/175. The 
plasma temperature T p is the same for all three cases. 
Notice that the electrons are strongly degenerate at all 
p and T displayed in Fig. At p > 4 x 10 10 g cm~ 3 
carbon nuclei cannot survive in dense matter because of 
beta captures; at p > 2 x 10 10 g cm -3 oxygen nuclei will 
also be destroyed by beta captures. Therefore, it is un- 
reasonable to extend the CO diagram to higher densities. 
We will be interested in nuclear fusion reactions 



need 



(A C ,Z C ), 



(9) 



where A c = A4 + Aj and Z c = Z^ + Zj refer to a com- 
pound nucleus c. For our example in Fig. [21 we have 
three reactions, C+C, C+O, and O+O ( Section |TT| . The 
experimental cross sections for these reactions show no 
very pronounced resonance structures and can be de- 
scribed in the framework of non-resonant reaction for- 
malism ^lUil as discussed above. Notice, that the data 
for the C+C and C+O reactions exhibit some resonant 
structures demonstrated in Fig. ^ I n the lower panel 
of this figure we show the Gamow-peak energy range 
as a function of temperature for these reactions in the 

We see 



thermonuclear burning regime (Section IIII Bfl 
that at high enough T the Gamow-peak range covers the 
energy range, where the oscillations are experimentally 
measured. However, at these values of T the Gamow- 
peak range is sufficiently wide and the energy integration 
in the reaction rate should smear out the oscillatory be- 
havior. For lower T the Gamow peak is narrower, but 
it shifts to low energies inaccessible to laboratory exper- 
iments. In the absence of experimental and theoretical 
data on the presence of oscillations at these low energies 
we will adopt the standard assumption that the reactions 
in question can be treated as non-resonant in applications 
to stellar burning. 

The shaded strips Fig. |2Jshow the T — p domains most 
important for these reactions. We will describe them in 
more detail in Section llVl 

To study a reaction © we introduce the so called ion- 
sphere quantities 



Z.Zje 2 



T (i) = ZjZje 



rp(l>) _ __ 



h ( AnZiZj 



fee V 2 ^i. 



1/2 



(10) 



where pij — m u AiAj / A c is the reduced mass of the re- 
acting nuclei, and n.y = 3/(47ra?-). Basing on the ion- 
sphere model of a strongly coupled Coulomb plasma, one 
expects (e.g., Ref. 0) that a,j characterizes an equi- 
librium distance between neighboring nuclei i and j, Fij 

describes their Coulomb coupling, is the tempera- 

(v) 

ture for the onset of strong coupling, and is a local 
Debye temperature (for oscillations of ions i and j). In 



an OCP we have a< 



a, Fi 



T, Ti, 



T p . We will also 



r-Bij = h 2 /(2p ij Z l Z j t 



(11) 



which reduces to the ion Bohr radius in the case of equal 
ions j = i. In addition, we will need the parameter 



A.i -+- A.j 



grgjj ( pX N {Z) ^ 1/3 
(Zl /3 + Z 1/3 ) V 2(A)m u 



AiAjZiZjizj /3 + Z X P) 
pX N (Z) 



1/3 



(12) 



_ (A) 1.3574 x 10 11 g cm" 3 , 

which corresponds to the parameter A introduced by 
Salpeter and Van Horn in the OCP case. 

In the following sections we will discuss nuclear burn- 
ing in MCP for the five burning regimes introduced in 
Ref. [a and analyzed in detail for OCP in our previous 
work [ll|. We will demonstrate that the formalism de- 
veloped for OCP can be adapted to more complex MCP 
scenarios. 



B. Classical thermonuclear reaction rate 



In the classical thermonuclear (weak screening) regime 
(T ^> Ty) the reacting ions constitute an almost ideal 
gas 0] • The rate for non-resonant fusion processes (such 
as considered here) is well known, 



R 



th 



l2Ef S{Ef) 



4 ni rij 
1 + Sij V 3pij knT 



exp(-Tjj): 



(13) 



where S(E) is the astrophysical factor; is the Kro- 
necker delta, which excludes double counting of the same 
collisions in reactions with identical nuclei (i = j); 
Efj = Tk-QTij/3 is the Gamow peak energy (the relative 
energy of colliding nuclei which gives the major contri- 
bution into the reaction rate) and 



Jl,. 72 ?2 A \ 1 / 3 



2k B Th 2 



(14) 



is the parameter which characterizes the Coulomb barrier 
penetrability. This parameter can be written as 

m = 3 (^/2) 2 / 3 (^/fc B T) 1 / 3 , E% = 2p lJ Z 2 Z 2 e A /h 2 . 

(15) 

Then 

R ^ = TTs- s{E ^ )r -¥ PthFt ^ (16) 

where reij is a convenient dimensional factor defined by 
Eq. (|11|) . i^th is the exponential function, and P t h is the 
pre-exponent, 

871- 1 / 3 / P a \ 2/3 
pH«). P lh = _(^) . (17) 
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Typically, the main contribution into the reaction rate 
comes from suprathermal particles (Ef^ ^> k^T), and 
the Coulomb barrier is very thick (r^- 3> 1). The reaction 
rate decreases exponentially with decreasing T. Typi- 
cal Gamow-peak energy ranges for the C+C, C+O, and 
O+O reactions in the thermonuclear regimes are shown 
in Fig. ^ These energies are defined as £J? k - AE < E < 



E 



P k 



AE, with AE 



EfHnT. 



C. Thermonuclear regime with strong screening 



The thermonuclear regime with strong plasma screen- 
ing operates in the temperature range T^- < T < 
where the plasma ions constitute a strongly coupled 
Coulomb system (liquid or solid). The majority of ions 
are confined in deep Coulomb potential wells but the 
main contribution into the reaction rate comes from a 
very small amount of highly energetic su pra thermal ions 
which are nearly free (see, e.g., Refs. 0,^3)- However, 
neighboring plasma ions strongly screen the Coulomb in- 
teraction between the reacting ions. The screening sim- 
plifies close approaches of the reacting ions, required for 
a Coulomb tunneling, and enhances thus the reaction 
rate (with respect to the classical thermonuclear reaction 
rate). 

In analogy with the OCP case (e.g., Ref. [13), the en- 
hancement can be included into the exponential function 

Ffh, 



F sc exp(-Ty), F sc = exp(hij), 



(18) 



where F sc is the enhancement factor and hij is a function 
of plasma parameters. 

We will analyze hij in the same manner as was done 
in Ref. for the OCP. For this purpose we notice 
that the reacting nuclei move in the potential W(r) — 
ZiZje 2 /r — Hij(r), where Hij(r) is the plasma potential 
created by neighboring plasma ions. In the thermonu- 
clear regime, (r) is almost constant along a Coulomb 
tunneling path. Accordingly, hij can be split into two 



terms, hi* = h 



- 1,(0) 



h 



(i) 



The main term iv/ is obtained 



assuming that Jiijir) ~ -ffij(O) is constant along a tun- 
neling path; a small correction is produced by a weak 
variation of the plasma potential along this path. We 
will discuss hf^ explicitly in this section and introduce 

/iy phenomenologically in Section IIII Gl when we pro- 
pose an analytic approximation for the reaction rate in 
all regimes. 

It is well known (e.g., Refs. 0) that h\f = 
Hij(0)/kQT, where fly(O) is the properly averaged 
(mean-field) plasma potential at r = 0. Thus defined, 
Hij (0) becomes a thermodynamic quantity which can be 
expressed through the difference of the classical Coulomb 
free energies for a given system and a system with the 
two reacting nuclei merging into a compound nucleus 



(e.g., Ref. 41]). However, a strongly coupled classical 
multi-component Coulomb liquid obeys very accurately 
the linear mixing rule (see Ref. H5 for recent results). 
Using this rule, one obtains 



fo(Ti) + f (Tj) - f (T c ), 



(19) 



where /o(r) is the Coulomb free energy per ion in an 
OCP (in units of k^T). This formula was derived by Jan- 
covici [42] for an OCP and generalized by Mochkovitch 
for an MCP (see Ref. 0). The function / (r) is very 
accurately determined by Monte Carlo simulations. For 
instance, according to DeWitt and Slattery |44| , the func- 
tion jfo(r) for a classical one-component Coulomb liquid 
at 1 < r < 170 can be approximated as 

/ (r) = -0.899172T + (1/s) 0.602249 V s 

-0.274823 lnT - 1.401915, (20) 

where s — 0.3230064. It gives a highly accurate expres- 
sion for , but it is inconvenient for an analytic in- 
terpolation of the reaction rate ( Section IIII Gfl . Instead, 
we will use a simpler linear expression /o(r) = — 0.9T 
provided by the ion-sphere model [14| . 



1,(0) s~isc p 



(21) 



C? 



0.9 



z 5/3 _ z $/3 _ z 5/3 



,1/3 



,1/3 



2^2 " j 



This expression seems crude but, actually, it is accurate. 
For a nuclear reaction in an OCP (Zi = Zj = Z c /2), 
this equation gives C sc = 1.0573, very close to the value 
1.0754 inferred 11] from Eq. ipDjl. In the range of 1 < 
r,j < 170 the ion-sphere model (|21|l gives the enhance- 
ment factor exp(ft, J -^' ) ) which is systematically lower than 
the more accurate enhancement factor, given by Eqs. I|19|) 
and l|20|l . The maximum difference of these enhancement 
factors for charge ratios 1/5 < Zi/Zj < 5 reaches w 15 
at the highest value of r^- = 170. This difference can 
be regarded as insignificant because at such the en- 
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74 



hancement factor itself is as huge as exp(h^) 
For lower the expression (|21|) is more accurate. For 
instance, for = 50 and Zi/Zj = 5 it underestimates 
the enhancement factor only by a factor of 3. In the range 
of 1 < Tij < 10 for 1/5 < Zi/Zj < 5 the underestimation 
does not exceed a factor of 1.5. Notice that in the MCP 
the factor CfJ depends on Zi and Zj. 

Although the above approach is more accurate, one 
usually calculates h\j by extrapolating the mean-field 
plasma potential Hij(r), obtained from classical Monte 
Carlo sampling, to r — > 0. In particular, Ogata et al. 
[45L 1461 ] used this method to study the enhancement 
of thermonuclear burning in the liquid phases of OCP 
and binary ion mixtures (BIMs). The leading term in 
Eqs. (19) and (20) of Ref. 0> equivalent to the leading 
terms (JT^J and is h[f = T l3 (1.148-0.00944 Inly- 
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0.000168 (In Tii ) 2 ). This expression was employed also 
by Kitamura |47| for constructing the analytic approx- 
imation for the reaction rates in OCP and BIMs in all 
reaction regimes (although by that time a more accurate 
expression was obtained by Ogata |48j for the OCP using 
path integral Monte Carlo simulations). Comparing the 
expression of Ref. |46l with (|19|) and l|21|) we see that the 
expression of Ref. :46] systematically overestimates the 
plasma enhancement. In the OCP the overestimation 
reaches a factor of ~ 40 for T = 170, which is not 
very significant. However, the coefficients in this expres- 
sion are independent of Zi and Zj , in disagreement with 
Eq. (|21|l . As a result, the overestimation increases with 
the growth of the charge ratio Zi/Zj, reaching ~ 150 and 
~2x 10 4 at Zi/Zj = 2 and Zi/Zj = 5, respectively, for 
Tij = 170. Therefore, when the difference of charges Zi 
and Zj increases, the results of Refs. and ^3 become 
less accurate. The nature of inaccuracy comes from the 
problems of extrapolation of Hij{r) to r — > 0. It was 
analyzed by Rosenfeld for the OCP 49] and the MCP 

Equations (|20[l and l|21|l become invalid in the regime 
of weak screening (T^ <C 1; Section llll Bp . where the well 
known Debye-Hiickel theory should be used. In partic- 
ular, the screening function in the weakly coupled MCP 
becomes [lil 



Z.Zje 1 
fc B Tr D 



srliz^zfz] 
(Z) 



1/2 



(22) 



where rrj is the ion Debye screening length and (Z 2 ) 
Y] j Z 2 Xj. Introducing fy, we obtain 



,(0) 



3/2 



Ci: 



'3(Z 2 )(zj /3 + Z) lz f 



1/2 



K -MJ-« . - ^ 8(Z)ZiZj j 

(23) 

For reactions in an OCP, we have C = \/3. In an MCP, 
Qj depends on ion charge numbers. 

A simple phenomenological interpolation which repro- 
duces the strong and weak screening limits (Eqs. I|21l) and 
(|23|I % 1 and combines them in the range from ^0.1 to 
~ 1 is 



4 0) = c-| 7 c r|/ 2 /[(c? 7 c /Ci J -) 4 + r|] 1 / 4 . 



(24) 



Because accurate calculations of the MCP free energy of 
ions in this range are absent, we cannot test the accu- 
racy of our interpolation. However, the plasma screening 
enhancement of reaction rates at these values of is 
weak and the interpolation uncertainty does not affect 
strongly the reaction rates. 



D. Zero-temperature pycnonuclear regime 



This regime operates at low temperatures, T < 
QST^ /\n{T$ at which thermal effects are neg- 



ligible and all the ions occupy ground states in their po- 
tential wells. The Coulomb barrier is penetrated owing 
to zero-point vibrations of ions around their equilibrium 
positions. Because the vibration amplitudes are gener- 
ally small, neighboring pairs of ions (closest neighbors) 
make the major contribution into the reaction rate. 

Generalizing Eq. (35) of Salpeter and Van Horn || to 
the MCP case we can present the pycnonuclear reaction 
rate as (see, e.g., Eq. (7) in Ref. [Ijj) 



1 + St. 



{ &ij Pij ) av 7 



(25) 



where i>ij is the number of nearest nuclei j around a nu- 
cleus i, is the reaction rate for a fixed pair ij, and the 
brackets (. . .) av denote statistical averaging over an en- 
semble of such pairs. For instance, it is currently thought 
that the OCP of ions at zero temperature forms a bec 
crystal. In this case, any ion is surrounded by 8 closest 
neighbors, with the equilibrium distance between them 
d = (37r 2 ) 1 / 6 a, where a is the ion-sphere radius. 

According to Eqs. (35), (37), and (39) of Ref. 0, the 
reaction rate for a pair of neighboring ions in an OCP 
can be written as 



P = D p 



A 3-C pl 5(£ P k) 



exp 



a, 



(26) 



where A and rs are given by Eqs. I|12l) and 111|) . respec- 
tively (for the OCP); while -D p k, C p i and C exp are con- 
stants which depend on a Coulomb barrier penetration 
model and on the lattice type. Finally, the characteristic 



reaction energy is E pk ~ hu> p — ksT p 



where uj p is the ion 



plasma frequency defined by Eq. JHJ (a typical frequency 
of zero-point ion vibrations). 

Salpeter and Van Horn Q used the WKB approxi- 
mation and considered two models of Coulomb barrier 
penetration in the bec lattice, the static and relaxed lat- 
tice ones, to account for the lattice response to a motion 
of tunneling nuclei. Later Schramm and Koonin 9] ex- 
tended this consideration taking into account the dynam- 
ical effect of motion of surrounding ions in response to 
the motion of the tunneling nuclei in the relaxed lattice. 
In addition, theyconsidered the fee Coulomb lattice. The 
results of Refs. HE1 are analyzed in Ref. 11]. Note that 
the pycnonuclear burning rates for fee or bec OCP crys- 
tals are very similar. We expect that for an amorphous 
OCP they are of the same order of magnitude. 

Pycnonuclear reactions in the MCP require compli- 
cated calculations (which, hopefully, will be done in the 
future). We will restrict ourselves to a simpler consid- 
eration based on similarity criteria and some general as- 
sumptions. Even a state of the MCP at low temperatures 
is not clear. It can be a regular lattice (with defects); a 
uniformly mixed state; an ensemble of phase-separated 
domains, etc. fSection IIII A|> . It can also be a combina- 
tion of these states. 
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1. Uniformly mixed MCP 

Let us start with a uniformly mixed MCP. An obvious 
generalization of Eq. 126|) to the MCP would be to replace 
r B ~ * r Bij an d A — » Ay in accordance with Eqs. (|llfl and 
(|12|l . This replacement should correctly reflect the rescal- 
ing of inter-ion distances and oscillator frequencies in an 
MCP within the ion-sphere model (see, e.g., Ref. 

However, the rescaling may be not exact. We will 
take this into account by adopting a simplified assump- 
tion that any pair of close neighbors behaves as an el- 
ementary oscillator with an equilibrium separation dij 
and an oscillator frequency wy. For an OCP we have 
A ~ h 2 /(Am n d 2 uj) 2 and lu ~ uj p . In the MCP we ex- 
pect to have A — ► Ay ~ ft 2 / '(2pijd 2 jUJij) 2 , where dij and 
ujij are the actual equilibrium distance and the effective 
oscillator frequency, respectively. With this replacement 
from Eq. I)2(:i[| we obtain 



Pij — ^pyk ' 



A 



S(Ef 



a 



hr hj 



exp 



cxp 



(Ay) 1 / 2 , 



(27) 



Now we assume that the actual values dij and can 
deviate from the ion-sphere values df^ and tof^ , and in- 



troduce the quantities 



Otdi 



dij/df* 



/ (0) 



(28) 



*'3/ "ij i —wrj ~ \l l ~%j 

which measure the deviations. They will be treated as 
parameters to be varied within reasonable limits. Then 

Ay = Ay a\ij, a Xi j = l/ia^al^), (29) 

where A y - is given by Eq. I|12fl . 

Another way to improve the ion-sphere rescaling for a 
pycnonuclear reaction in a BIM was proposed by Ichi- 
maru et al. |5l| who used the formal relation Ay ~ 
rBij /d^ and suggested that Ay = Ay/ady, which cor- 
responds to a\ij — l/a>dij- Therefore, they allowed dij 
to be different from dfy but assumed that the oscillator 
frequency wy adjusts to this new separation following 
the ion-sphere rescaling, so that a 2 ^- = l/ajjy. In con- 
trast, we allow the separations and oscillator frequencies 
to deviate independently, and our consideration is more 
general. 

Following Ogata at al. (see their Eq. (7)), we as- 
sume that the number of closest neighbors j around the 
ion i in Eq. I|25|l is (i^y)av = Snj/n (which is appropri- 
ate for a uniform mix). The factor '8' may be approxi- 
mate but it affects only the pre-exponent of the reaction 
rate which is much less significant than the exponentially 
small probability of Coulomb tunneling. Substituting 
(1271) and J2HJ| into (gSJ, we obtain 



R pyc 



D,. 



Ariifij 



i(Z) 



1 + % (z 1/3 + z 



.^ L - exp 

A,"? 1 



a 



cxp 



(Ay) 1 / 2 



(30) 



This equation has the same structure as Eq. l|16f) and 
can be written as 



R pyc 



UiUj 

1 + Si, 



S(Efj ) j ^pyc-Fpyc- 



(31) 



with 



P, 



i(Z) 4D pyc 



Co 



pyc -( Z ^ + zfr a^ ' pyc 

For numerical evaluations, we have 



(Ay) 1 / 2 



(32) 



RT = io 46 c pyc 



8pX N XiXjAiAj (A)Z 2 Z 2 
(1+^yW 




cm" 3 s~\ (33) 



where C pyc = D pyc /(8 x 11.515); the density p is ex- 
pressed in g cm~ 3 and the astrophysical factor S(Ef.) 
is in MeV barn. The reaction energy is Ef- ~ hoja — 



(p) 



y —*3 

, . The main parameter regulating the reaction 

For sufficiently 



rate is Ay in the exponent argument 
low densities, Ay is very large, strongly suppressing the 
Coulomb tunneling. With growing p the barrier becomes 
more transparent and the reaction rate increases. 

For the OCP with ot\ = 1 Eqs. (gU||-@ reduce to the 
well known equations for zero-temperature pycnonuclear 
burning in a crystalline lattice. The constants C pyc , C p i 
and C cxp , obtained using various techniques, have been 
analyzed in Ref. • In Table [H] we present these pa- 
rameters for three models. The first model is optimal 
(seems to be the most reliable). It is the static-lattice 
model of Salpeter and Van Horn 8] for the bcc crystal. 
The second and third models are phenomenological; they 
have been proposed in Ref. [ijj , basing on the results of 
Refs. 0, 0- 03- D3 ■ The second model gives the upper 
limit of the reaction rate and the third gives the lower 
limit (for both - bcc and fee - crystals). We expect that 
the reaction rate in an amorphous OCP would lie within 
the same limits. 

Returning to an MCP, we must additionally specify 
the scaling factor a.\y in Eq. (|2"§|> . Some attempts have 
been made |45ll5l| to determine proper inter- ion separa- 
tions dij (i.e., the values of ayy) at T = in BIM solids 
from positions of first correlation peaks in radial pair dis- 
tribution functions of ions, gij{r). These functions have 
been calculated by Monte Carlo sampling. Such studies 
require powerful computer resources and may be incon- 
clusive at present. This is clearly seen from similar (and 
simpler) attempts of the same group |45ll4ril5lll52l | to de- 
termine deviations of inter-ion separations from the ion- 
sphere scaling in BIM liquids (using first-peak positions 
of gij(r)). The authors applied their results to study 
an effect of these deviations on thermonuclear burning in 
the strong screening regime. Our analysis of those results 
shows that no statistically significant deviations from the 
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ion-sphere scaling have been found (and no associated ef- 
fects on nuclear burning can actually be predicted) . This 
conclusion is strengthened by the critical analysis of these 
works by Rosenfeld [491 IBoj ] . Therefore, no reliable infor- 
mation has been obtained on the violation of the ion- 
sphere scaling of inter-ion separations in the BIM liquids 
and solids (on the level of a few percent, which is most 
likely the real uncertainty in the determination of gij (r) 
peak positions in the cited publications). 

Equally, proper oscillator frequencies Wjj (and param- 
eters a^ij) in BIMs could be determined from molecular 
dynamics simulations but such simulations have not yet 
been performed. 

In the absence of precise microscopic calculations of 
dij and Wy we naturally assume that the optimal values 



are a di 



a\ij 



1. In order to maximize the 



reaction rate we propose to increase axij (somewhat ar- 
bitrarily) by 5%, and in order to minimize the rate we 
propose to reduce it by 5%. These proposed values are 
also listed in Table ITU Notice that it is difficult to expect 
that the variations of adij and d U y in Eq. I|28|) are fully 
independent. An increase in the inter- ion separation dij 
should cause a decrease in the oscillator frequency w^ ; 
these variations should be partially compensated in the 
factor a M j = V^L^tf)- 



2. Regular MCP lattice 

Now we turn to the case a regular MCP lattice, which 
can be drastically different from a uniform mix. The 
central point is the availability of closest neighbors ij. 
If they are absent, the reaction (0 occurs via Coulomb 
tunneling of more distant ij pairs and becomes strongly 
suppressed. 

The closest-neighbor condition depends on the crystal 
type. For instance, consider a binary bcc crystal com- 
posed of ions i and j. There are eight pairs of closest 
neighbors in a basic cubic cell (formed by one ion in the 
center of the cell and any other ion in a vertex) . If all ions 
are of the same type (e.g., Xi — 1), then all eight pairs 
participate in the same reaction ii (the OCP case). In 



the BIM case (x. 



4) we have an ion i in the cen- 



ter of the cell and ions j in vertices, and all eight pairs of 
closest neighbors participate only in the reaction ij. The 
reactions ii and jj will be strongly blocked because any 
pairs ii and jj are not closest neighbors. Then the ion i 
in the center of the basic cell will be able to react with 
six ions i in centers of adjacent cells. The equilibrium 
distance between these pairs is a factor of 2/V3 « 1.155 
larger, than between the closest neighbors, which will ex- 
ponentially suppress the ii reaction rate. 

One can construct more complicated MCP lattice 
structures and formulate appropriate blocking condi- 
tions. The strongest blocking of pycnonuclear reactions 
is expected in a regular crystal with many components. 
In that case the probability to find specified closed neigh- 
bors ij could be very selective. 



3. Other MCP structures 

If the matter consists of domains of separated phases, 
pycnonuclear burning occurs mainly within these do- 
mains. The rate of the reaction ii in domains containing 
ions i can be calculated from Eq. I|33(l . assuming OCP 
(xi = 1), and then diluted by a volume fraction occu- 
pied by the phase i. The burning of different ions i and 
j occurs on interfaces between corresponding domains. 
Roughly, in a BIM with spherical bubbles containing N 
less abundant ions (i or j) the OCP reaction rate ii can 
be diluted with respect to the rate in a uniformly mixed 
state by a factor of 1 /TV 2 / 3 (because the reactions occur 
on interfaces). 

Pycnonuclear burning in a lattice can be drastically af- 
fected by lattice impurities and imperfections (see, e.g., 
Ref. 0). We expect that the effect of impurities or im- 
perfections can be included in Eq. (|33|) if we treat them as 
members of the MCP. Because the burning is mainly reg- 
ulated by the parameter , the rate should be extremely 
sensitive to variations of equilibrium distances and oscil- 
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FIG. 3: (color online) Rates of pycnonuclear C+C, C+O and 
O+O reactions versus density for the optimal burning model. 
The dot-and-dashed lines show C+C burning in a pure carbon 
crystal and O+O burning in a pure oxygen crystal. Other 
lines are for the CO mixture (xc = 0.5); the solid lines refer 
to a uniform BIM; the long-dashed lines are for a regular CO 
bcc crystal; the dots are for phase-separated matter. The 
inset shows the same curves for the C+C reaction on a larger 
scale; the thicker short-dashed lines give the maximum and 
minimum reaction rates for the uniform CO mixture in the 
ion-sphere model (ax = 1); the thin short-dashed lines are 
the same but allowing for the variation of ax (see the text for 
details) . 
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lation frequencies for the reacting nuclei (see Eq. 
Small variations can induce exponentially huge jumps or 
drops of the Coulomb tunneling probability (for instance, 
in response to the decrease or increase of inter-ion sepa- 
rations). The impurities and imperfections may be rare 
but give the leading contribution into the reaction rate. 

Moreover, different pairs ij in an MCP may be ex- 
posed to different local conditions and have different sep- 
arations and oscillator frequencies. One can incorporate 
these effects by introducing the averaging (. . .) av over an 
ensemble of such pairs in Eq. i|25|) . 

For illustration, in Fig.[3|we show the density depen- 
dence of C+C, C+O and O+O pycnonuclear reactions 
in carbon-oxygen BIMs. The astrophysical factors are 
taken from Section [n] In the main part of the figure 
we use the optimal reaction model from Table [H] The 
dot-and-dashed lines show the C+C burning in a pure 
carbon matter (xq — 1) and O+O burning in a pure 
oxygen matter (xc — 0) while other lines are for the CO 
mixture with xq = \- The solid lines present the re- 
action rates in the uniform CO mixture. For the C+C 
and O+O reactions these lines go slightly lower than the 
dash-and-dot lines because of the reduced amount of car- 
bon and oxygen in a BIM (as compared to pure carbon 
or oxygen matter). The long-dashed lines give the re- 
action rates in the bcc CO regular crystal. The C+O 
burning in this crystal has the same rate as in the CO 
uniform mix (the solid line), in our model. In order to 
illustrate the blocking effect (produced by the absence of 
closest CC or 00 neighbors in the crystal) we have cal- 
culated the C+C and O+O reaction rates using Eq. 
but increasing the distances between the reacting nuclei 
in Eq. Ij28(l by a factor of a ah — 2/V3, where i=G or O 
(see above). Because the local oscillator frequencies for 
these CC and 00 pairs are unknown, we have assumed 
the ion-sphere rescaling (a^ u — l/a\ u ). The blocking 
effect is strong, reducing the reaction rates by 5-7 orders 
of magnitude. Nevertheless, even with this blocking, the 
C+C reaction is much faster than the C+O and O+O 
ones because of much lower Coulomb barrier. Finally, 
we have also calculated the BIM reaction rates assum- 
ing phase separation of C and O. We do not plot the 
respective C+C and O+O reaction rates because they 
are almost indistinguishable from the corresponding dot- 
and-dashed lines (the rates are only twice lower than in 
the respective OCPs). However, the C+O burning in this 
case can be strongly suppressed because it occurs only at 
interfaces between separated phases. For example, the 
dotted line shows such a C+O burning rate assuming 
separation into domains which contain 10 6 ions. 

The inset in Fig. [3] displays the C+C reaction rate on 
a larger scale. The solid and dot-and-dashed lines are 
the same as in the main part of the figure. Short dashed 
lines present the maximum and minimum reaction rates 
for the uniform CO mixture; they reflect uncertainties 
of existing theoretical models for pycnonuclear burning. 
These lines are calculated using the parameters from Ta- 
ble [n] Thicker short-dashed lines give the maximum and 



minimum reaction rates using the ion-sphere model (ne- 
glecting deviations from the ion-sphere rescaling; i.e., as- 
suming ciacc = !)• Thinner short-dashed lines are the 
estimated maximum and minimum reaction rates tak- 
ing into account possible deviations from the ion-sphere 
model. The current theoretical uncertainties of the reac- 
tion rate are really large. 



E. Thermally enhanced pycnonuclear regime 



This regime operates in the temperature range 
QJoT$/\d.(T$/T$) < T < 0.5 T$ (see Ref. i). At 
these temperatures, the majority of nuclei occupy ground 
states in their potential wells but the main contribution 
into the reaction rate comes from a small amount of nu- 
clei which occupy excited bound states. This regime 
is difficult for theoretical studies. We will follow the 
approach of Ref. ^l| > which is based on an analytical 
approximation of the WKB calculations performed by 
Salpeter and Van Horn [8|. We generalize this approach 

"by 



to an MCP in the same manner as in Section fill Dl 
rescaling A and T p in accordance with Eqs. I|12fl and ( HUP . 
In this case 



^ yc ( T ) 

< c (o) 



-i = 



7TJ2 



exp 



T- 

-A y 



T 



-AT*' /T 



(34) 

...„, , , ].; - a u ijT-?', while SI, fix, and A 
are model-dependent dimensionless constants. We adopt 



where = hujij/k^ 



0.95 to maximize the reaction rate and a u 



1.05 to minimize it fTable lHjl . In analogy with Ref. |ll| 
the characteristic energy of the reacting nuclei can be 
taken in the form 



E 



P k 



ZiZjt 



exp 



-A 



(35) 



The first term is the reaction energy in the zero- 
temperature pycnonuclear regime, while the second term 
describes an increase of Ef* with growing temperature 
in the thermally enhanced pycnonuclear regime. Equa- 
tion (|35|l is approximate (based on the results of Ref. @ 
as explained in |ll| ) but we expect that it correctly re- 
flects the main features of the accurate expression (to 
be derived in future, more elaborated calculations). The 
analogous expression (26) for OCP in Ref. 0] contains 
two new free factors, C\ and C2, in the first and sec- 
ond terms, respectively. They were set C\ = C2 = 1 in 
numerical calculations 0] because the theory Q is insuf- 
ficiently precise to determine them. We do not introduce 
these factors here to avoid additional unknowns, but they 
could be introduced in the future. An uncertainty in Ef^ 
should not greatly affect the reaction rates. 

After Salpeter and Van Horn @ the thermally en- 
hanced pycnonuclear burning in an OCP was studied by 
Kitamura and Ichimaru |53| assuming that the reacting 
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TABLE II: Coefficients in the interpolation expressions for a reaction rate for the optimal model of nuclear burning and for the 
models which maximize and minimize the rate. The parameters Ct, a uij ar e different for MCP and OCP (the values for 

OCP nj are given in brackets). For an MCP, the models assume a uniformly mixed state (see the text for details). 



Model 


Cexp 


Cpyc 


C p i 


Ct 






A 


Optimal 


2.638 


3.90 


1.25 


0.724 (0.724) 


1 (1) 


i (i) 


0.5 


Maximum rate 


2.450 


50 


1.25 


0.840 (0.904) 


1.05 (1) 


0.95 (1) 


0.35 


Minimum rate 


2.650 


0.5 


1.25 


0.768 (0.711) 


0.95 (1) 


1.05 (1) 


0.65 



nuclei move in an angle-averaged, radial static mean-field 
potential determined from Monte Carlo sampling of clas- 
sical Coulomb systems. Although this approach is less 
justified than the WKB approximation of Ref. 0, the 
results are in a reasonable agreement (see Ref. for de- 
tails). Kitamura 01 generalized the results of Ref. |5^| 
to the case of BIMs using the ion-sphere rescaling rule 
{a\i 3 = a^ij = 1). 

The thermally enhanced pycnonuclear burning is as 
sensitive to the microphysical structure of the MCP 
as the zero-temperature pycnonuclear burning (Section 
1111 D|) . The above comments refer to a uniformly mixed 
MCP. In the case of a regular MCP lattice the reaction 
will be suppressed by the same blocking effects as dis- 
cussed in Section IlIIDI The case of phase separation has 
the same features as at T = 0. Illustrative examples will 
be given in Section 



F. The intermediate thermo-pycnonuclear regime 



This regime is realized at temperatures TJy /2 < T < 

TJy which separate the domains of classical and quantum 
motion of the reacting nuclei. The calculation of the 
reaction rate in this regime is most complicated. We 
will describe this rate by a phenomenological expression 
presented below. 

The reaction is mainly determined by the nuclei which 
become slightly unbound and can move freely through 
the dense matter, fusing not only with the closest neigh- 
bors (pycnonuclear regime), but with other nuclei (ther- 
monuclear regime). We expect that the transition from 
the pycnonuclear to the thermonuclear regime with the 
growth of temperature in a uniformly mixed MCP is suf- 
ficiently smooth. When the number of freely reacting nu- 
clei becomes large, the dependence of the reaction rate on 
the details of the MCP microstructure should disappear. 



G. Single analytical approximation in all regimes 

Our phenomenological expression for the temperature 
and density dependent reaction rate, which combines all 
the five burning regimes and assumes a uniformly mixed 
MCP at low temperatures, is a straitforward generaliza- 



tion of the expression for the OCP considered in Ref. [Tl| . 

R lJ (p,T) = m c (p) + AR t3 (p,T), 



AR t] ( Pl T) 



rijUj S(Ef-) 



l + Sij h 



T, 



F = exp -nj + C sc fy ^e" AT - /T - A -g- , 



T 



p= 8^/3 (E ?^ 



V3 2V3 y k T 



(36) 



In this case, v = yffi /[((%/& ) 4 + T?// 4 ; R™ c {p) 
is the density dependent pycnonuclear reaction rate at 
zero temperature discussed in Section IIII Dl ARij (p, T) 
is the density and temperature dependent part including 
an exponential function F and a pre-exponent P. The 
quantities ry and Fy are similar to the familiar quantities 
Tij and Tij , but contain a "renormalized" temperature T, 



(! 



2/3 / E: 



v = 



Z l Z j e 2 

~ 1 

aijksT 



k B T 

f = 



1/3 



r 2 + c'2(r l ( f ) ) 2 , (37) 



where Ct is a dimensionless renormalization parameter 
specified below. Equations l|36l) and l|37|l are analogous 
to Eqs. (27) and (28) of Ref. \u\. The term ARij (p, T) 

and the renormalized temperature T are introduced to 
match the equations in thermonuclear and pycnonuclear 
regimes. The renormalized temperature reflects the fact 
that the thermal energy k^T of plasma ions in the ther- 
monuclear case is replaced by a temperature- independent 
zero-point energy in the pycnonuclear case. 

(p) ~ 
For high temperatures T 3> T£_- we have Tij — > Ty, 

~ —> T. In this case AR lj {p,T) -> 
» R*j"(p), and Eq. reproduces the 

thermonuclear reaction rate (Sections IIII Bl and IIII C|l . 
At low temperatures T < the quantities r,;, 

the quantum temperature" T$ 



r 



H&(j>,T) » R pyc 



ij 



and T contain "the quantum temperature" 
determined by zero-point ion vibrations, rather than 
the real temperature T. In the limit of T — > 
we obtain f lj = l/[(Xij ) 1/2 (727T) 1 / 6 C T ] and T tj = 

3(7T/Ay) 1/2 /(2 7/6 4 /3 )- 
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Following Ref. [ll| we require that at T < T^f the 
factor exp(— Ty ) in the exponential function F, Eq. 1)36(1 . 

reduces to exp(— C OX p/(Ai :) ) 1 / 2 ). This would allow us to 
obey Eq. l(51|l by satisfying the equality 



3^/(2 7 / 6 4 /3 Hc cxp K 2J r 1/2 - 



(38) 



The double exponent factor in F, Eq. (|36[1 . will corre- 
spond to the double exponent factor in Eq. I|34|l . Tak- 
ing C exp and a\ij from Table ITT1 we can determine Ct- 
These parameters are also listed in Table [H] In the MCP 
they are different from those in the OCP because in the 
MCP we introduce an additional parameter a\ij (Section 
Imp)) . The values of C T for the OCP [HI are given in 
Table ITT1 in brackets. 

Finally, the quantity 7 in Eq. (|36|l and the reaction 
energy Ef* in the astrophysical factor S(Eff) can be 
chosen in the same way as in Ref. [TTj] . 



(r 2 7l 



(7f) 2 7 2 ) / (T 2 + {f^f 



(39) 



EfY = Kb. 



(p) 



AT, 



(p)- 



exp 



T 



(40) 

where 71 = 2/3 and 72 = (2/3) (C p i + 0.5). 

Thus, we propose to use the analytic expression i|36|) 
for the reaction rate in a uniformly mixed MCP with 
the following parameters: 

(1) The parameter C sc of strongly screened thermonu- 
clear burning is given by Eq. (|21() . 

(2) The parameters C cxp , C pyc , C p \, ot\ij of zero- 
temperature pycnonuclear burning, and the parameters 
a^ij, A and Ct of thermally enhanced pycnonuclear 
burning are given in Table ITT1 

In this way we obtain (Table [HJ three models for any 
given non-resonant nuclear fusion reaction @ in a uni- 
formly mixed MCP. One is the optimal model, the second 
gives the maximum reaction rate, and the third the min- 
imum reaction rate. For the OCP, it is sufficient to set 
= oiujij — 1 1 which reduces the present results to the 
results of Ref. [llj • 

The uncertainties of the reaction rate become larger if 
a cold MCP forms a regular lattice or undergoes a phase 
separation or contains impurities and defects. All these 
cases can be approximately taken into account in the 
same way as discussed in Section IIII Dl For instance, a 
reaction in a regular lattice can be strongly blocked by 
the absence of closest reacting neighbors (Sections IIII Dl 
and IIII El also see Section iTvTl . 

Our formula for a uniformly mixed MCP gives a 
smooth behavior of the reaction rate as a function of 
temperature and density, without jumps at the solidifi- 
cation point (in analogy with an OCP, see Ref. [H). In 
the cases of other MCP microstructures such jumps may 
appear. 



Our formula is flexible. Its parameters could be tuned 
when new microscopic calculations of reaction rates ap- 
pear in the future. Moreover, the formula can be im- 
proved even if a new information on MCP properties (not 
on reaction rates directly) appear in the literature (for in- 
stance, on the deviations from the ion-sphere scaling at 
T = 0). 

More complicated expressions for the reaction rates in 
the OCP and uniform BIMs were proposed by Kitamura 
|47| . His ex pres sions are mainly based on the results 
of Refs. HI lH H |H i| (in the different regimes) 
which are not free of approximations (see Ref. |ll| for 
details). His expressions for BIMs are obtained assum- 



ing the ion-sphere rescaling rule (a\i 



1) and 



are, therefore, more restricted than our expression. Their 
derivation implies that they are valid for uniformly mixed 
BIMs. In particular, they do not take into account block- 
ing effects in regular binary lattices. 

In contrast to our formula, Kitamura took into account 
the effects of electron screening (finite polarizability of 
the electron gas). However, these effects are relatively 
weak; their strict inclusion in the pycnonuclear regime 
is complicated. We do not include them but, instead, 
take into account theoretical uncertainties of the reaction 
rates without electron screening. The results of Kitamura 
[47| for an OCP lie well within these uncertainties. His 
results for BIMs in the thermonuclear regime with strong 
screening and essentially different charges of reacting nu- 
clei are less accurate than our results (Section IIII C|) . 



IV. NUCLEAR BURNING IN A 
CARBON-OXYGEN MIXTURE 

For illustration of our results, we analyze nuclear re- 
actions in a dense 12 C 16 mixture. Figure 01 shows 
the temperature dependence of the C+C reaction rate 
at p = 5 x 10 9 g cm -3 (for the same reaction models 
as in Fig. |3J). The dot-and-dashed line is the optimal 
model (Table [Hj for a pure carbon matter (from Ref. 
|ll|'). Other lines are for CO mixtures with equal num- 
bers of C and O nuclei (xc = \)- The solid line is the 
optimal model for the uniform mixture; the thicker short- 
dashed lines show the maximum and minimum reaction 
rates in such a mixture in the ion-sphere approximation 
(«acc = o^cc — 1); the thinner short-dashed lines are 
the same but beyond the ion-sphere approximation (note 
that variations of a^cc appear to be much less impor- 
tant than variations of a\cc)- The long-dashed line is for 
a CO regular lattice (in the same approximation which 
has been used in Fig.|3J). A sharp jump of this curve is 
associated with the melting of the crystal (Fig. 0) , which 
destroys the blocking of the C+C burning and amplifies 
the reaction rate. 

Figure0]shows the reaction rates in all burning regimes 
(except for the classical thermonuclear burning which 
would require higher temperatures; see Fig. |2J. The 
horizontal parts of the curves for logT[K] < 7.7 refer 
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log 10 T [K] log p I0 [g cm- 3 ] 



FIG. 4: (color online) The C+C reaction rate versus temper- 
ature at p = 5 x 10 9 g cm -3 . The dot-and-dashed line is the 
optimal model for a pure carbon matter; other lines are for 
a CO mixture with xc = 0.5. The solid line is the optimal 
model for a uniform mixture; the long-dashed line is for the 
regular bcc CO crystal at low temperatures (the vertical part 
indicates the melting point). Thicker short-dashed lines give 
the maximum and minimum reaction rates for the uniform 
CO mixture in the ion-sphere model (a\ = 1); thinner short 
dashed-lines are the same but allowing for deviations from the 
ion-sphere model (see the text for details). 



FIG. 5: (color online) Carbon ignition curves in 12 C 16 mat- 
ter. The dot-and-dashed line is the optimal model for carbon 
burning in pure carbon matter. The solid and dotted lines are 
optimal models for uniform CO mixtures with xc —0.5 and 
0.1, respectively. Other lines are for CO BIMs with xc = 0.5. 
The short-dashed lines give the highest and the lowest theo- 
retical ignition curves for uniform mixtures. The long-dashed 
line is for the CO bcc crystal at low temperatures. 



to the zero-temperature pycnonuclcar burning; the re- 
spective reaction rates are independent of T as discussed 
in Section IIII Dl and displayed in Fig. |3| The temper- 
ature range 7.7 < logT[K] < 8.3 corresponds to the 
thermally enhanced pycnonuclear regime. The reaction 
rate starts to grow up with increasing T (Section IIII Ej l. 
The rate remains highly uncertain because of the same 
reasons as in the zero-temperature pycnonuclear regime. 
The next temperature range 8.3 < logT[K] < 8.6 cor- 
responds to the intermediate thermo-pycnonuclear burn- 
ing (Section llll El) . Theoretical uncertainties of the reac- 
tion rate become smaller. Finally the temperature range 
logT[K] > 8.6 refers to the thermonuclear burning with 
strong plasma screening. The theoretical uncertainties 
become much smaller although the enhancement of the 
reaction rate by the plasma screening effects is huge; with 
increasing T this enhancement weakens and the reaction 
rate matches the classical thermonuclear rate (see Fig. 6 
of Ref. .11]). The presence of oxygen slightly reduces the 
C+C reaction rate (by reducing the amount of carbon 
nuclei at a given density). 

Based on our expression for the reaction rates, in Fig. [21 
we plot the T — p domains (shaded strips), where the 
C+C, C+O and O+O reactions are most important. The 



domains for the C+C and O+O reactions are presented 
for pure carbon or pure oxygen matter [xc — 1 and 
xq = 0, respectively). For the C+O reaction we have 
taken the CO mixture with xc = \- A domain for any 
reaction ij is restricted by two lines, along which the 
characteristic burning time Tj = rii/Rij of nuclei i is con- 
stant (taken to be 10 6 years for a lower line and 1 year 
for an upper line, for example). Above the upper line 
the reaction ij is so fast, that the nuclei i cannot sur- 
vive for a long time. Below the lower line the reaction is 
so slow that the nuclei i survive almost forever. There- 
fore, the strips represent the temperature and density 
domains of greatest relevance for the carbon and oxygen 
nucleosynthesis through the reactions under discussion. 
For determining these domains, we have taken the opti- 
mal model from Table [HI The domains do not change 
significantly under variations of fractional numbers of C 
and O within reasonable limits. For densities p < 10 9 
g cm -3 the strips are almost horizontal; nuclear burn- 
ing in them proceeds in the thermonuclear regime and 
the reaction rates depend mainly on the temperature. In 
contrast, the strips become almost vertical at low temper- 
atures reflecting the pycnonuclear burning regime where 
the reaction rates depend mainly on the density. 

The strips show a strong heterogeneity of the different 
reactions. It is evidently caused by different heights of 
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Coulomb barriers. With increasing p and/or T in the 
CO matter, carbon will burn first in the C+C reaction 
and could be burnt almost completely before reaching the 
T — p domain, where the C+O reaction can be efficient. 

Finally, we have studied the carbon ignition curve, 
which is a necessary ingredient for modeling nuclear ex- 
plosions of massive white dwarfs (supernova la events) 
and carbon explosions in accreting neutron stars (super- 
bursts). The ignition curve is usually determined as the 
line in the T ' — p plane (Fig. EJ) , where the nuclear energy 
generation rate equals the local neutrino energy losses. 
At higher T and p (above the curve) the nuclear energy 
generation exceeds the neutrino losses (which cool the 
matter) and carbon ignites. We have calculated such 
curves for CO mixtures. All reactions (C+C, C+O, and 
O+O) have been taken into account but the C+O and 
0+0 reactions have appeared to be unimportant owing 
to the heterogeneity of nuclear burning. The presence of 
oxygen affects carbon ignition only through the C+C re- 
action rate and the neutrino emission rate. The neutrino 
energy losses have been assumed to be produced by plas- 
mon decay and by electron-nucleus bremsstrahlung. The 
neutrino emissivity owing to plasmon decay has been ob- 
tained from extended tables calculated by M. E. Gusakov 
(unpublished); they are in good agreement with the re- 
sults of Itoh et al. 54]. The neutrino bremsstrahlung 
emissivity has been calculated using the formalism of 
Kaminker et al. [5j|, which takes into account electron 
band structure effects in crystalline matter. For a CO 
mixture, this neutrino emissivity has been determined 
using the linear mixture rule. 

The dot-and-dashed line in Fig.[S]shows the carbon ig- 
nition curve, calculated using the optimal model of car- 
bon burning in a pure carbon matter (xc — 1; from Ref. 
[ll|). The solid and dotted lines are the same curves in 
CO mixtures with xc = 0.5 and 0.1, respectively (assum- 
ing the optimal reaction model and a uniform mixture at 
low temperatures). At p < 10 9 g cm -3 the curves de- 
pend weakly on the density because carbon burns in the 
thermonuclear regime. At T < 10 8 K the curves depend 
weakly on the temperature because carbon burns in the 
pycnonuclear regime. A strong bending of the curves in 
the density range from ~ 10 9 g cm~ 3 to ~ 3 x 10 9 g cm~ 3 
is associated with the transition from the thermonuclear 
burning to the pycnonuclear one. With decreasing the 
carbon fraction, the ignition curve shifts to higher T and 
p, mainly because of the decrease of the C+C reaction 
rate. 

The short-dashed lines in Figure [S] show the uncer- 
tainty of the solid ignition curve (xq = \ , a uniform CO 
mixture) associated with the uncertainties of the reac- 
tion rates (assuming the maximum and minimum reac- 
tion rates from Table [HJ . In the thermonuclear regime 
the uncertainties are small, while in the pycnonuclear 
regime they are substantial. The long-dashed line shows 
the ignition curve calculated under the assumption that 
a regular bec CO lattice is formed in the CO mixture 
(xc = 5) after the crystallization. The blocking of the 



C+C reaction rate by oxygen ions in the bec lattice shifts 
the ignition curve to higher p. The sudden break of this 
line is associated with the crystallization (analogous to 
the break of the long-dashed line in Fig. 0}. 

The carbon ignition curve obtained by equating the nu- 
clear energy generation and the neutrino losses becomes 
unreliable for T < 10 s K (e.g., Ref. [H). The main rea- 
son is that this curve falls in the T — p domain, where 
the characteristic carbon burning time is unrealistically 
large (exceeds the Universe age). In addition, the neu- 
trino emission becomes very slow, inefficient for carry- 
ing away the nuclear energy; thermal conduction can be 
much more efficient. As a result, the carbon ignition con- 
dition becomes nonlocal, complicated, and dependent on 
a specific model (a neutron star or a white dwarf, etc.). 



V. CONCLUSIONS 

We have studied the problem of Coulomb barrier pen- 
etration for non-resonant nuclear fusion reactions in a 
dense MCP of atomic nuclei. We have considered all five 
nuclear burning regimes (Sections IIII BllrTlF|l and ana- 
lyzed calculations of nuclear reaction rates in an MCP 
for these regimes, available in the literature. We have 
proposed (Section IlII G|) a unified phenomenological ex- 
pression for the reaction rate valid for all regimes. It gen- 
eralizes an analogous expression proposed recently |ll| 
for an OCP. The expression contains several parameters 
which can be varied to account for current theoretical 
uncertainties of the reaction rates. 

Our main conclusions are: 

1. The reaction rates in the thermonuclear regimes 
(with weak and strong plasma screening) can be 
calculated sufficiently accurately. In the regime of 
strong screening and for reacting nuclei with non- 
equal charges, our expression is more accurate than 
those proposed in the literature (Section [ill C|) . 

2. The reaction rates in other regimes (zero- 
temperature and thermally enhanced pycnonu- 
clear regimes; intermediate thermo-pyeno nuclear 
regime) are much less certain. They are very sen- 
sitive to currently unknown microphysical correla- 
tion properties in an MCP (a uniform mix, a regu- 
lar crystalline lattice, a phase separated matter, a 
matter with impurities and defects); they are much 
richer in physics than in the OCP case. 

3. At low temperatures, we have mainly considered re- 
actions in a uniform mix. Other MCP microstruc- 
tures can strongly decrease or increase the reaction 
rates. For instance, the reactions in a regular MCP 
lattice can be strongly suppressed by the absence 
of nearby reacting nuclei. 

4. Our phenomenological expression can be improved 
(Section IIII G|) after new calculations of the reac- 
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tion rates or main properties of the MCP are per- 
formed. It would be important to know the ac- 
tual microstructure of the MCP at low tempera- 
tures (first of all, the availability of closest neigh- 
bors, local separations and oscillation frequencies 
of neighboring nuclei, particularly in the presence 
of impurities and lattice defects). 

5. Although our main formula in Section IIII CI as- 
sumes a uniform mix at low T, the presented results 
are sufficient to understand qualitatively the reac- 
tion rates for other cases (following prescriptions of 
Section llTTD)l . 

For illustration, we have considered ('Section llVfl C+C, 
C+O, and O+O nuclear reactions in a dense carbon- 
oxygen mixture, that is important for the structure and 
evolution of massive white dwarfs (supernova la explo- 
sions) and accreting neutron stars (as sources of super- 
bursts) . For this purpose we have calculated and parame- 
terized the appropriate astrophysical factors ( Sect ion ITT1). 
The main results of our analysis are as follows: 

1. The ranges of densities and temperatures, where 
C+C, C+O, and 0+0 reactions are most impor- 
tant, look like narrow regions in the temperature- 
density diagram (Fig. 0); the regions do not 
strongly overlap which means strong heterogeneity 
of these reactions. 

2. With increasing density and/or temperature, car- 
bon starts burning first in the C+C reaction (be- 
cause carbon nuclei have lower Coulomb barrier); 
this reaction is most important for the nuclear evo- 
lution of CO mixtures. 



3. Carbon burning in the C+C reaction is affected 
by the presence of oxygen. The effect is simple in 
the thermonuclear regimes but more complicated 
in other regimes (at low temperatures). 

4. Carbon ignition in a CO mixture occurs (Fig. [SJ in 
thermonuclear regimes as long as p < 10 g cm~ 3 
(and T > 3 x 10 8 K). It can be calculated quite 
accurately. With decreasing the fraction of carbon, 
the ignition curve shifts to higher p and T. 

5. At p > 10 9 g cm" 3 and T < 3 x 10 8 K the ignition 
condition becomes uncertain fSection llVfl . The for- 
mation of a regular CO lattice after the crystal- 
lization can block the C+C reaction and shift the 
carbon ignition to higher densities. 

Our consideration in this paper was general. More 
quantitative nuclear network simulations involving ther- 
monuclear and pycnonuclear burning in dense stellar 
matter are currently in progress. 
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